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INTRODUCTION 


It  ha?  long  been  known  that  the  ionosphere  at  high  latitudes  often 
displays  density  irregularities  and  structured  electric  fields  which  are 
of  various  magnitudes  in  scale  and  strength  [Fejer  and  Kelley,  1980].  The 
recent  resurgence  in  interest  in  high  latitude  irregularities  has  been  due 
to  the  observation  of  scintillations  by  DNA’s  Wideband  Satellite  in 
association  with  the  gradient  in  the  total  electron  content  (TEC)  at  these 
latitudes  [Fremouw  et  al.,  1977;  Rino  et  al . ,  1978].  A  variety  of 


experimental  techniques  have  now  been  employed  to  study  the  structure  at 
high  latitudes  and  its  morphology  is  beginning  to  be  understood  better. 
On  the  basis  of  incoherent  scatter  radar  returns,  the  plasma  density 
contours  have  revealed  the  existence  of  high-latitude  enhancements  which 
have  scale  sizes  transverse  to  and  along  the  magnetic  field  on  the  order 
O'  hundreds  of  kms  and  have  frequently  been  referred  to  as  "blobs" 
[Vickrey  et  al.,  1980;  Kelley  et  al.,  1982].  Simultaneous  Incoherent 
scatter  and  VHF  scintillation  measurements  were  used  by  Vickrey  et  al. 
(1980)  to  establish  the  co- location  of  -  1  km  scale  irregularities  and  the 


F  layer  blobs.  Ground-based  optical  and  digital  ionosonde  measure-ments, 
in  conjunction  with  the  particle  measurements  on  the  DE-2  satellite,  have 
found  the  presence  of  large  scale  plasma  irregular lties  (called  "patches") 
at  polar  cap  latitudes  [Weber  et  al . ,  1984],  Various  mechanisms  have  been 
proposed  for  the  cause  of  these  large-scale  plasma  irregular  it ies ,  such  as 
generation  by  particle  precipitation  and  convection  from  other  regions 
[Kelley  et  al . ,  1982;  Kelley  and  Vickrey,  1984;  Weber  et  al.,  1  984,  Sojka 
and  Schunk,  1986].  Associated  with  the  large  scale  blobs,  there  are 
density  gradients  at  "walls"  of  these  blobc  (primarily  in  the  horizontal 
east-west  direction)  that  are  observed  to  be  structured  [Vickrey  et  al . , 
1980;  Rino  et  al.,  1983,  Basu  et  al.,  1984;  Cerisier  et  al.,  1985,  Tsunoda 
et  al.,  1985].  It  is  generally  believed  that  the  structure  at  the  so- 


called  "intermediate”  scale  sizes  (-  1-10  kms)  is  responsible  'or  th* 


observed  satellite  scintil latior.s  . 
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A  piasma  blob  at  high  latitudes  Is  reminiscent  of  artificially 
released  clouds  In  the  Ionosphere  and  this  .'act  has  motivated  us  to  do 
analogous  analysis  of  the  L  «  _3  Instability  Tor  the  blob  situation.  In 
this  connection,  we  mention  that  the  £  «  3  instability  has  been  applied  to 


the  naturally  occurring  transverse  density  gradients  at  mid-latitudes 


[Scar.napieoo  et  al . ,  1975]  and  to  the  high-latitude  situation  [Xeskir.en 
and  Ossakow,  1982;  Keskinen  and  Ossakow,  1983;  Mitchell  et  al.,  1985].  In. 


addition,  we  note  here  that  yet  another  gradient-driven  instability,  the 
current  convective  Instability  (which  occurs  in  presence  of  a  parallel 
current)  has  been  invoked  to  interpret  high-latitude  structure  in 
situations  which  may  be  _E  «  B_  stable  [Ossakow  and  Chaturvedi,  1979; 
Chaturvedi  and  Ossakow,  1979;  Huba  and  Ossakow,  1980;  Chaturvedi  and 
Ossakow,  1981].  Most  of  the  past  work  for  the  E_  *  B_  and  current- 
convective  instability  was  carried  out  in  the  approximation  of  an 
infinitely  long  blob  parallel  to  the  magnetic  field.  These  earlier 
studies  clearly  need  to  be  modi'ied  for  the  blob-related  structure,  since, 
as  has  been  noted  earlier,  the  blobs  have  finite  length  in  the  direction 


parallel  to  the  field-lines.  The  work  on  the  _E  x  Instability  in  barium 
clouds  has  shown  that  the  growth  rate  is  severely  modified  for  the  long 
wavelength  modes  when  the  finite  parallel  size  of  the  cloud  is  considered 
[Volk  and  Haerendel,  1971;  Goldman  et  al.,  1976;  Sperling  et  al.,  1989; 
Drake  et  al.,  1985;  Drake  and  Huba,  1986].  A  similar  analysis  carried  out 


for  the  current  convective  instability  has  demonstrated  that  the 


instability  growth  rate  is  reduced  on  the  order  of  two  orders  of  magnitude 
for  the  scintillation  causing  scale  sizes  (for  the  observed  blob  sizes) 
[Huba  and  Chaturvedi,  1986]. 


In  this  paper  we  investigate  the  _E  *  B  Instability  'or  the  high 
latitude  blob-like  situation.  Thus,  in  our  analysis -we  have  included  the 
effects  of  the  finite  parallel  blob  size,  and  ion-inertia  (which  can 


become  important  at  high  altitudes  where  the  ion-neutral  collision 
frequency  is  small  [Ossakow  et  al.,  1°78]).  We  'ind  that  the  growth  rate 
for  the  E  *  3  instability  is  modi'ied  by  the  finite  parallel  size  of  the 


blob,  with  the  longer  wavelengths  showing  more  severe  reductions  in  growth 


rate  than  the  small  and  intermediate  scale  sizes.  The  effect  of  ion- 


inertia  at  high  altitudes  also  leads  to  a  reduction  in  the  growth  rate, 


out  the  reduction  is  relatively  moderate  for  the  small  and  intermediate 
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wavelengths  U  <  10  km)  but  Is  substantial  at  the  long  wavelength?  (A  >  10 
km).  Further,  the  mode  1?  localized  in  the  blob-region  parallel  to  the 
magnetic  field.  This  has  some  notable  consequences.  For  Instance,  the 
mode  may  "“main  relatively  unaffected  by  the  regions  afar.  Thus,  the 
blob-associated  structure  of  F-region  altitudes  may  develop  even  in  the 
presence  of  a  conducting  E- region.  The  mapping  of  structure-associated 
potential  'ields  is  also  modified  due  to  the  localization.  The  electro¬ 
static  potential  is  roughly  constant  within  the  enhancement  (cloud)  but  as 
it  maps  into  the  background  ionosphere,  the  mapping  distance  is  finite  due 

to  the  dissipation.  This  mapping,  originally  discussed  by  Farley  [1959], 

1/2 

follows  the  relationship  (L  /L  ]  *  (<j  /<j  )  ,  where  L  is  the  character- 

II  II 

lstic  field  scalesize,  a  is  the  conductivity  of  the  medium,  and  the 
subscripts  |,  refer  to  the  parallel  and  transverse  components  to  the 
magnetic  field,  respectively.  We  find  that  this  mapping  may  not  extend  to 
large  parallel  distances  for  structured  fields  at  higher  altitudes 
where  v^n  (and  thus  o  )  is  very  small.  We  discuss  these  results  in  the 
final  section. 

The  plan  of  the  paper  is  as  follows.  In  the  next  section,  we  present 
the  basic  equations  and  briefly  discuss'  the  equilibrium.  In  the  third 
section,  the  nonlocal  mode  equation  is  derived.  The  fourth  section 
contains  the  numerical  results.  Finally,  in  the  fifth  section  the  results 
are  discussed. 


II.  BASIC  EQUATIONS  AND  EQUILIBRIUM 

A  set  of  nonlinear  equations  describing  the  evolution  of  a  plasma 
density  enhancement  has  been  derived  in  Sperling  et  al.  [1984]  and  Drake 
et  al.  [1985].  We  briefly  present  the  derivation  here.  A  plasma  blob  is 
assumed  to  have  finite  length  along  the  uniform  magnetic  field  B_  =*  Bg 
z  with  a  background  neutral  wind  -  Vn  x.  The  plasma  is  assumed  to  be 
cold  (i.e.,  T  -  0);  see  Appendix  A  for  a  discussion  of  finite  temperature 
effects.  The  basi..  equations  are 
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In  the  above  n  and  v  are  the  density  and  'laid  velocity,  respectively, 
a  -a 

o'  species  a,  v  is  the  a  species  -  neutral  collision  frequency,  v  is 
an  - 

the  electron-ion  collision  frequency,  is  the  ion-electron  collision 
frequency  and  a  refers  to  electrons  (e)  or  ions  (i).  The  electric  and 
magnetic  fields  are  represented  in  terms  of  potentials  as 

.  3A  _ 
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—  0  z 


where  $  is  the  electrostatic  potential  and  A2  is  the  vector  potential 

associated  with  the  magnetic  field  produced  by  the  self-consistent  plasma 

currents.  We  consider  only  A_  since  J  >>  J  and  assume  that 
-  z  I  1 

,7Az  x  z\  «  BQ. 

The  electron  motion  transverse  and  parallel  to  Bq,  is  described  as 


v  «  -  -  7  d>  x  z 
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where  \>a  -  +  von  and  V  »  b  •  7,  b  •  3/|Bq|.  Also,  with  vo/flo  <<  1 

and  v,  <  n. ,  f  *  m  v  /m.v.  <<  1.  The  ion  dynamics  is  given  by 

in  -  i  e  en  i  in 
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where  -  (1  ♦  Vin/2I)  '  and  vin  “  vin  +  ^t.  In  01),  both  Pedersen 
and  Hall  responses  to  the  electric  field  and  neutral  wind  have  been 
included.  Mote  also  that  we  have  included  the  ion  inertia  Q/3t  in  v^n) 
in  the  above,  even  though  the  full  convective  derivative  has  not  been 
included  (v  •  7v  term  is  not  taken  into  account).  The  reason  for  doing  so 
is  the  fact  that  this  term  only  introduces  a  real  frequency  to  the 
solution  and  does  not  affect  the  stability  criterion. 

We  present  the  equations  derived  by  Sperling  et  al.  [1984]  for 
the  E  *  B  instability  with  the  ion- inertial  effects  included.  These 
equations,  derived  from  (3)— (12)  and  (1),  (4)  and  (5),  are 
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where  ne  -  mov0/ne  is  the  parallel  resistivity  and 
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Equation  (13)  ia  the  electron  continuity  equation,  (14)  is  the  charge 
conservation  equation  and  (15)  is  the  Ampere's  law.  Now  we  shall  discuss 
the  equilibrium  of  these  equations  briefly. 
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The  equilibrium  of  a  cold  plasma  enhancement  ha?  been  discussed  by 
Sperling  et  al.  [1984],  A  two-dimensional  cloud  localized  along  and 
across  the  magnetic  field  _Bq,  =*  n^(x,z)  with  *  0  for 
I x |  i  Xq  and  \z\  i.  zQ  is  described  by  equations 
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In  the  above,  it  has  been  assumed  that  the  background  ionosphere  is 
uniform  throughout  the  region,  and  xQ  and  zQ  are,  respectively,  the 
typical  transverse  and  longitudinal  scales  of  the  plasma  enhancements. 
From  (19)  —  (21  )  we  find  that  the  equilibrium  is  given  by  Az  *  <j>  -  0  with 
n(x,z)  an  arbitrary  function.  For  simplicity  we  consider  the  plasma 
density  to  be  given  by 


n(x,z)  -  n1 
n(x,z)  -  n^(x) 
n(x,z)  -  n_ 


z  <  -z. 
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III.  LINEAR  ANALYSIS  AND  DISPERSION  RELATION 

We  perturb  the  set  of  equation?  (13)~(15)  over  the  equilibrium 
discussed  in  the  previous  section  and  linearize  the  equations  with  respect 
to  perturbations.  The  perturbed  quantities  are  assumed  to  vary  a? 
f  -  f(z)  exp(Yt  +  k  y).  The  above  set  then  leads  to 
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The  above  equations  are  the  same  as  those  discussed  in  the  work  of 
Sperling  et  al.  [1984]  with  the  exception  of  the  second  term  in  paren¬ 
thesis  on  the  left  hand  side  of  (24)  (y/v,  )  which  reoresents  the  effects 

in 

of  ion  inertia.  The  following  definitions  have  been  used  in  (23)  and 


(24);  Yn  »  V  /L  ,  L  *  (3  In  n,/3x)  ! ,  D  -  v  .  c2/w  2,  u>  ^  =*  4irnr.e<Vm  , 

0  n  n  n  v  i  r  ei  pe  pe  0  e 

a  »  Q  Q./v  ./v.  ,  and  v.  /Q.  <<  1  was  assumed.  We  discuss  some  well-known 
e  i  el  in  m  l 

results  from  (23)  and  (24)  in  the  local  analysis  in  Appendix  A.  Now  we 
present  a  derivation  of  the  nonlocal  dispersion  relation. 

The  method  of  solving  (23)  and  (24)  we  follow  here  is  similar  to  that 
used  in  Sperling  et  al.  [1  934],  Drake  et  al.  [1985],  and  Huba  and 
Chaturvedi  [1986].  In  the  outer  regions  (outside  the  plasma  blob),  the 
ambient  ionosphere  is  assumed  to  have  no  density  gradient  and  a  decaying 
solution  is  assumed  for  the  modes  in  these  regions.  We  shall  distinguish 
between  the  two  outside  regions  by  their  differing  ambient  parameters 
(density,  v.  ,  etc.).  We  refer  to  the  region  below  the  blob  as  the  region 
1  (bottomside)  and  the  one  above  it  as  the  region  2  (topside).  In  these 
regions,  the  solution  is  written  as 
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The  inverse  of  decay  lengths  of  the  mode-associated  “ields  in  the  ambient 
ionosphere  are  given  by 

2  (l  +  — ) 

P-1.2  (26) 


where 
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In  the  above  p  -  1,2  refer  to  the  two  outer  regions  as  noted  earlier,  the 

term  (y/v.  j  in  (26)  describes  the  effects  of  ion  inertia  and  R  is  a 
in 
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measure  of  the  electromagnet! c  e'fects.  For  the  electrostatic 
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(k  D  >>  YJ,  R  »  1  and  for  the  electromagnetic  ease  i Y  >>  k  S  y ,  H  * 
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Inside  the  blob,  the  mode  is  represented  by  a  plane  wave  sol  at  i or. , 


Ik,  z 


i  f<  .  2 


A  .e 
2  i 


-Z  <  2  <  2„, 

J  VJ 


with  the  parallel  wavenumber  of  the  mode  inside  the  blob  being  given  cy 


Y(l 

7 


<  CQ  ' 


where  the  quantities  are  evaluated  in  the  blob  region  and  the  ior.-irert iai 
effects  are  represented  by  the  third  term  in  the  numerator  of  (29).  A 
solution  is  constructed  from  (29)  for  the  mode,  using  the  two  roots  of 
(k,  corresponding  to  +  roots), 

L*  ~ 

~  ik,+2  .  ik  2 

A.  -  A  e  1  ♦  B  e  (30 

A  similar  solution  is  written  out  for  The  solutions  are  matched  at 

the  boundaries  of  the  inner  and  outer  regions  to  derive  the  dispersion, 
relation.  The  appropriate  matching  conditions  are  obtained  by  integrating 
(23)  and  (2*0  across  the  boundaries  [z  -  +  2Q);  it  is  ‘‘ound  that  i  and 
A^  are  continuous  at  z  •  *  zQ.  Using  (25)  and  (30)  in  (23)  and  (2k)  with 
these  matching  conditions,  we  find  the  dispersion  equation,  is 

ii  + 

-i  k<  +  ki 

Vo"?13"  --V.77  *  J*  ;3,! 

'  - 

k, 

where  m  la  an  integer.  The  largest  growth  rate  occurs  for  m  =  0  node  and 
we  consider  only  this  mode  henceforth.  The  wavenumbers  k.  ^  and  k..  are 

t  ,  d. 

defined  in  (28)  and  (29),  nespecti veiy.  Equation  (31)  is  solved  numeri¬ 
cally,  but  first  we  shall  discuss  some  simple  analytical  results  ‘'non  it. 


:f  1 3cal  analysis  car.  be  obtained  from  (31)  in  the  limit 
.nere  '  »  yri  ,  ir.  which  ease  one  may  write  tan  (2k,  z  )  * 
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Tne  'i'.ite  length 

if  the  oi ib  and 

the 

a  ^e  : -.ct .  ir  if  tne  E  «  3  ir.staciiity  growth  rate.  We  have  noted  earlier 
that  we  distinguish  tne  ior.ospherio  regions  above  and  below  the  blob  by 
tneir  differing  ambient  parameters.  We  shall  also  distinguish  three 
separate  case-3  m  whien  the  effects  of  ion  inertia  are  introduced 
dif 'erectly  see  Tables  *  and  2).  In  case  I,  both  the  topside  and  the 
cottomsi  ie  ionosphere,  as  well  as  the  blob  correspond  to  ionospheric 
altitudes  where  tne  ion-inertial  effects  are  neglected.  In  the  case  II, 
the  ion  inertial  effects  are  Introduced  only  in  the  topside  ionosphere 
region  and  in  tne  case  III,  both  the  blob  and  the  topside  ionosphere  are 
considered  at  altitudes  where  ion-inertial  effects  may  become  important. 
We  present  some  simple  analytical  expressions  for  the  growth  rate  of 
the  E  <  3  instability  'or  the  above  three  cases  in  Appendix  B.  It  is 
shown  there  that  the  growth  rate  is  reduced  for  a  finite  length  blob  and 
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:he  previous  section.  The  parameters  used  have  been 


ijz  ,  z2  -  a,L  2,  and  L2  -  D  /yn. 
3  r  r  i  r  i  r  rO 


. '  u  ,  ~t  —  'j,  ,  u  ,  ,  '.i  i ..  i  *  i  he  blob 

y  r  3 

has  beer,  assumed  to  be  the  same  In  the  three  cases  ard  is 
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in  the  linear  local  theory  growth  rate  Y^-  10~^  sec-1 


taken  to  be  300  kms.  Similarly,  the  gradient  scale  length,  LN,  associated 
with  the  "walla"  of  blob  la  taken  to  be  uo  kms  and  the  neutral  wind 
velocity,  Vn,  la  taken  to  be  400  m/aec  in  all  the  three  cases,  resulting 

Most  of  the  otner 

parameters  and  symbols  are  explained  in  Tables  1  and  2  and  will  oe  stated 
when  introduced  below.  In  the  following  the  growth  rate  expressions  have 
been  computed  as  a  function  of  wavenumber,  the  blob  length  and  the  ion 
collision  frequency  for  various  cases.  In  addition,  we  have  also  plotted 
the  wave  function  within  and  outside  the  blob  for  various  oases,  to 
display  the  localization  of  the  mode  around  the  instability  region  (i.e., 
the  blob).  For  reference  we  note  that  1  km  modes  correspond  to  <^Lr  - 
14.5  and  10  km  modes  correspond  to  kyLr  -  1.45. 

Figure  1  shows  the  normalized  growth  rate,  Y/y  ,  plotted  as  a 

function  of  normalized  transverse  wavenumber,  k  L  ,  for  the  three  cases 

’  y  r 

mentioned  above.  In  case  I  (curve  A),  the  bottoms ide,  blob  and  topside 
are  all  taken  to  be  collisional  and  ion  inertial  effects  are  unimportant. 
For  large  ky  (short  transverse  wavelengths),  the  growth  rate  can  be  seen  to 

approach  asymptotically  the  local  theory  growth  rate,  Y  -  Y_ .  For  small  kv 

o  y 

(long  transverse  wavelengths),  the  growth  rate  is  reduced.  The  inter¬ 
mediate  scale  size  regime  corresponding  to  the  1  -  10  km  is  of  special 
interest,  for  these  irregularities  are  associated  with  the  signal 
scintillations.  It  can  be  seen  that  the  finite  blob  length  has  little 
effect  on  the  1  km  scalesize  mode,  where  growth  rate  is  close  to  the  local 
theory  value  (Y  “  0.9  YQ)  but  at  the  10  km  scalesize-er.d,  the  growth  rate 
is  down  by  almost  an  order  of  magnitude  (Y  *  0.2  Y^ ) .  When  inertial 
effects  are  introduced  for  the  topside  (case  II,  curve  3),  the  growth  rate 
curve  shows  an  overall  decrease,  and  the  reduction  at  long  wavelengths 
becomes  more  pronounced.  Thus  the  1  km-end  of  intermediate  scale  size 
shows  small  reduction  in  growth  rate  (Y  *  .75  y_),  while  the  '0  km  wave- 
length  modes  suffer  severe  reduction  (Y  *  .05  Y„)  in  growth  rate.  This 
trend  continues  for  the  case  when  the  inertial  effects  are  important  for 
both  the  blob  and  the  topside  (case  III,  curve  C).  Now  ever,  the  local 
growth  rate  (in  the  inertial  limit)  is  different  from  the  value  in  the 

collisional  case  (Y  *  0.6  Y_  in  the  inertial  case,  while  y  -  y.  in  the 

U  0 
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oollislonal  case).  The  '  xm-scales  sho w  less  reduction  in  the  growth  rate 
due  to  the  “inlte  Clod  length  eT'ects  (y  *  0.35  y  )  than  the  ’0  km-size 
scale?  (Y  *  0.025  Y  ^ ) ,  a?  expected. 

In  Figs.  2a  -  2b  the  explicit  eT“ect  of  “ir.ite  blob  size  is  ahowr.  or. 
the  growth  rates  of  1  <m-  and  iG  xm-sl-e  modes  Tor  the  cases  I  and  III, 
respectively.  Figure  2a  refers  to  case  I,  and  we  see  that  the  ’  km  scale 
sizes  (cjrve  A)  experience  lesser  reductions  in  the  growth  rate  than  the  '0 
xm-size  modes  (curve  3;.  It  can  be  seen  that  'or  very  long  blobs 

( approaching  tne  infinitely  long  blob  limit,  or  the  local  theory  limit), 
the  growth  rate  Tor  k^L^  “  "4.5  (i  -  1  Km)  approaches  the  local  theory 
value  (Y  «  7,)  and  'or  small  blob-lengths,  the  growth  rate  is  progressively 
reduced.  This  pattern  is  more  pronounced  Tor  the  longer  wavelengths 
( 'K,j-r  *  1  -10  kms )  Tor  which  the  growth  rate  ir.  the  in.Tinite  blob 

limit  (Y  »  0.3  yq)  is  less  than,  the  local  theory  value,  Y  »  YQ  (which  is 

valid  only  'or  short  wavelengths  such  that  kyLN  >>  1).  A  similar  trend  is 

'ound  'or  growth  rate  behavior  as  a  “unction  oT  blob-length,  when  the 
inertial  e““ects  are  Introduced.  This  is  shown  in  Fig.  2b,  and  the 
inertial  e““ects  have  been  included  Tor  both  the  blob  and  the  topside 
regions.  The  growth  rates  Tor  both  the  scale  sizes  U  -  1,  10  *m)  are 

reduced  in  this  (inertial)  case  "or  the  inTir.ite  olob  limit  “root  the  value 
correspond]  n.g  to  the  ooLlisional  limit.  And  as  the  blob  lengths  get 
smaller,  the  reduction  in  tne  growth  rate  is  more,  with  longer  wavelengths 
(curve  3)  showing  more  pronounced  reduction  than  the  small  wavelengths 
( curve  A) . 


In  Figs.  3  -  5  we  have  plotted  the  wave  “unction?  Tor  the  modes  within 
the  blob  (the  instability  region)  and  outside  it  (the  ambient  ionosphere) 
to  show  the  localization  oT  the  mode?  in  the  biob  region  and  its  extension 
into  the  ambient  ionosphere.  In  Fig.  3.  the  wave  “unction  is  plotted  “:r 
the  short  wavelength  mode  (k^G^  *  ‘b.5,  \  -  i  xm).  Figures  3a  and  3b 
refer,  respectively,  to  cases  I  and  II.  It  can  be  easily  seen  that  in  the 
collisior.al  bottomside  region  the  mode  extends  to  smaller  distances  ir.  the 
ambient  ionosphere  thar.  the  topside  ionosphere  where  collisions  are  low. 
This  mapping,  dependent  on  the  ambient  coll  is ior.al  ity ,  shows  “or  tne 
collisior.al  topside  (Fig.  3a)  the  mode  extending  to  ^mailer  distances  than 


The  mapping  of  modes  into  the  ambient  ionosphere  i?  also  scale  sice 

dependent.  This  is  shown  in  the  Figs.  4  and  5.  In  Fig.  4,  wavenumbers 

k  L,  “  ’4  5  (Fig.  4a),  and  k  L  ■  1.45  (Fig.  4b)  are  plotted  for  the  case 
y  r  y  r 

when  all  the  regions  are  taken  to  be  collislor.al  (case  I,  Table  2). 
Clearly,  the  longer  wavelengths  extend  'arther  into  the  background  than  the 
smaller  wavelengths.  In  Fig.  5,  which  corresponds  to  the  case  in  which  the 

topside  is  assumed  inertial  (ease  III,  Table  2),  this  pattern  is  repeated, 

with  the  potential  'ields  mapping  to  larger  distances  on  the  topside  (which 
is  less  collisioral)  than  on  the  bottomside  (which  is  more  coll  is  ior.al ) . 
Again,  as  one  would  expect,  this  mapping  is  to  larger  distances  for  the 
longer  wavelength  modes  (Fig.  5b,  k^L^  *  1.45)  than  for  the  small 

wavelength  mode  (Fig.  5a,  k^L^  ■  14.5). 

An  interesting  point  to  note  from  the  Figs.  3  -  5  is  that  in  the 

highly  collisional  bottomside  region,  the  smaller  scales  do  not  map  to 
large  distances.  Thus  we  see  that  the  1  km-size  mode  maps  to  less  than  the 
blob  length,  zQ/2  (-  75  kms )  in  the  collisional  bottomside  (Figs.  3a  and 
4a).  This  implies  that  this  mode,  depending  on  the  altitude  of  the  blob, 
may  not  be  affected  by  a  conducting  E-region.  A  conducting  E-region  is, 
however,  likely  to  affect  the  development  of  longer  wavelengths  (>  10  kms), 
as  is  evident  from  Figs.  3b  -  4b,  where  the  modes  extend  to  several  times 
the  olob  length  (>  400  Kms). 

V.  DISCUSSION 

We  have  studied  the  effects  o'  finite  blob-length  on  tne  interchange 
instability  applicable  to  high  latitudes  and  have  included  the  e'fects  of 
'ir.ite  lor.  inertia.  The  experimental  observations  indicate  that  the  high- 
latitude  blobs  are  of  "inite  extent  parallel  to  the  magnetic  'ield  with 
typical  scalesizes  of  a  few  hundreds  of  km  [Tsunoda  et  al . ,  1985].  This 
means  that  they  span  regions  of  varying  plasma  parameters  in  altitude, 
e.g.,  the  collision  'requencies,  density,  etc.  At  high  altitudes,  the  ion- 
neutral  collision  frequency,  v,n>  decreases  and  therefore  the  ion-inertial 
e"eots  need  to  be  included  for  the  interchange  instability  analysis 
(OssaKow  et  al . ,  1978].  The  introduction  of  'inite  parallel  blob-size  in 


the  analysis  means  that  one  should  consider  the  '‘inite  parallel  wavelength 
effects  (k  »  o)  on  the  development  of  the  instability  and  it  further 
implies  that  the  coupling  of  the  blob  region  to  the  ambient  ionospheric 


regions  should  be  taken  into  account. 


There  are  some  obvious  similarities  of  dynamics  between  an  artificial 
plasma  cloud  and  a  plasma  blob  in  the  ionosphere.  3oth  represent  a  plasma 
density  enhancement  above  the  background,  both  are  of  finite  dimensions  and 
in  both  the  cases  their  motion  with  respect  to  the  background  results  in  a 
steepened  backside  or  a  wall  gradient  which  is  observed  to  be  striated  or 
structured.  Considering  these  common  features,  we  have  applied  the  theory 
of  interchange  instability  developed  for  plasma  clouds  by  Sperling  et  al. 
[1984]  and  Drake  et  al.  [1985]  to  the  blob  situation,  and  have  included  the 
ion-inertial  effects  in  the  analysis.  It  may  be  pointed  out  that  there  are 
also  some  differences  between  the  two  situations.  The  clouds  are 
relatively  small  in  parallel  size  (tens  of  kms)  as  compared  with  the  blob 
sizes  (a  few  hundred  kms)  and  have  a  sharper  transverse  density  gradient. 
The  constituent  ions  of  clouds  are  also  different  (usually  barium)  from  the 
background  ions  while  in  the  case  of  the  blobs,  the  two  are  the  same.  Due 
to  the  large  sizes,  blobs  occupy  regions  of  diverse  ionospheric  properties 
in  altitude.  Thus  the  plasma  properties  at  the  bottom  and  the  top  of  a 
blob  may  be  very  different  in  terms  of  collisionality,  temperature, 
density,  etc.  This  also  has  an  important  implication  on  the  profile  of  the 
blob  (as  opposed  to  that  of  a  cloud  in  whose  case  this  size  is  quite  small, 
and  thus  the  properties  of  the  two  ionospheric  regions  adjacent  to  the 
cloud  are  not  as  different).  Even  though  we  have  considered  a  'top  hat' 
model  of  the  blob  in  our  analysis  for  simplicity,  the  real  blob  would  be 
asymmetrical  in  altitude,  with  Its  boundary  in  the  ccllisional  bottomside 
remaining  sharp  over  the  time-scale  of  the  structure  development  as  opposed 
to  the  topside  where  this  boundary  would  be  relatively  more  diffuse  with 
respect  to  the  background. 


The  analytical  and  numerical  work  on  the  stability  of  blobs  with 
respect  to  the  interchange  mode  presented  in  the  previous  sections 
demonstrate  that,  as  in  the  case  of  the  plasma  clouds,  the  finite  parallel 
blob-ler.gth  results  in  an  overall  reduction  of  the  growth  rate  of  the 
Interchange  instability  from  the  Infinite  blob  case  value.  This  reduction 
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is  proportional  to  the  wavelength  of  the  mode:  the  longer  the  wavelength, 
the  larger  the  reduction  in  its  growth  rate.  In  addition,  the  effects  of 
ion-inertia  (important  at  high  ionospheric  altitudes)  are  found  to  reduce 
the  interchange  instability  growth  rate  somewhat.  The  consideration  of 

finite  parallel  blob-length  implied  that  the  effects  of  finite  k^  be  con¬ 
sidered  on  the  interchange  instability.  The  parallel  extent  of  the  modes 
is  proportional  to  the  transverse  wavelength  in  ionosphere  and  this  rela- 
tionshlp  is  roughly  given  by  the  relation,  (L^/L  )  *  ^a|/ojJ'  ’  wiiere  L  13 
the  characteristic  scale  size,  o  is  the  ionospheric  conductivity,  and,  j_ 
and  |  represent  components  perpendicular  and  parallel  to  the  magnetic  field 

respectively.  Clearly,  as  one  goes  higher  in  altitude,  o  (which  is  pro- 

X  -i 

portional  to  v  )  decreases  and  a  (which  is  proportional  to  vg^) 

increases.  Thus,  with  increasing  altitude,  a  given  transverse  scale  size 

mode  extends  to  larger  and  larger  parallel  distances.  This  fact  explains 

the  more  pronounced  reduction  of  interchange  instability  at  longer 

transverse  wavelengths  for  the  finite-length  blob  situation.  As  the  mode 

extends  out  of  the  olob  (the  instability  region)  into  the  ambient 

ionosphere,  it  attempts  to  move  plasma  electrodynamically  in  these  regions 

(where  the  density  is  uni"orm  and  this  represents  a  drag  for  the  system). 

As  a  result,  the  growth  rate  of  the  interchange  process  is  reduced.  For 

longer  transverse  wavelengths,  the  parallel  extent  of  the  mode  is  also 

larger,  and  thus  the  'load*  also  higher  for  the  circuit,  and,  therefore, 

the  growth  rate  reduction  is  more.  Another  consequence  of  the  relationship 
1/2 

(L  /L  ]  «  (o  /a  is  the  fact  that  the  smaller  and  intermediate  scale- 
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sizes  are  localized  to  shorter  parallel  distances  around  the  blob  (the 
instability  region).  Hence,  depending  on  the  altitude  where  the  blob  is 
located,  this  localization  may  be  small  enough  so  that  the  mode  does  not 
extend  to  the  E-region  on  the  bottoms ide.  This  would  mean  that  the  E- 
region  would  have  little  impact  on  the  development  or  the  structuring 
process  in  the  high-latitude  F-region,  whether  or  not  the  E-region  is 
conducting.  The  satellite  signal  scintillations  are  generally  believed  to 
be  caused  by  the  intermediate  scale  sizes  (few  100  m  -  few  km)  and  for  a 
set  of  typical  parameters,  we  have  found  in  our  numerical  results  that  this 
may  indeed  be  the  case  [see  Figs.  3a,  3b>  ^a  and  5a]. 


The  long  parallel  extension  of  modes  on  the  topside  also  means  that 
the  regions  of  low  neutral  density  (and  therefore  low  v,R)  are  being 
traversed  by  these  modes.  Thus  ion-inertial  effects  need  to  be  included 
and  this  has  an  effect  of  reducing  t:.e  interchange  instability  growth  rate 
[Ossakow  et  al . ,  1978].  This  reduction  is  also  wavelength  dependent,  the 
larger  transverse  wavelengths  experience  higher  reductions  in  the  inter¬ 
change  growth  rates  (see  Figs.  1  and  2,  curves  B  and  C,  and  also  the 
nonlocal  growth  rate  expressions  (35)  and  ( B6 ) ) .  However,  the  interchange 
process  may  proceed  rather  unimpeded  at  short  and  intermediate  scales,  with 
the  reductions  in  the  growth  rates  less  than  an  order  of  magnitude  even  at 
higher  (-  450  kras)  altitudes,  a  result  in  agreement  with  Ossakow  et  al. 
[1973]  (see  Figs.  1  and  2). 

Thus  we  find  that  the  interchange  process  associated  with  the  high 
latitude  'inite-length  blobs  is  severely  hampered  fir  long  wavelengths 
(>  10  km)  due  to  finite  blob-size,  possible  coupling  to  a  conducting  E- 
region  on  tne  bottomside,  and  due  to  inertial  effects  as  modes  couple  to 


higher  altitudes  on  the  topside.  On  the  other  hand,  the  small  ar.d 
intermediate  scale  sizes  (<  1  km  and  1-10  km,  respectively)  may  still 
experience  the  process  of  Interchange  unhindered  for  the  observed  blob 
lengths  (-  300  km)  for  the  F-region  altitudes  (i.e.,  at  altitudes  around 
400  -  500  km).  There  is  evidence  that  the  h.gh-latitude  blobs  may  be 
displaying  such  structures  (3ino  et  al . ,  1983;  Cerisier  et  al . ,  *985; 
Tsunoda  et  al . ,  1985).  Therefore  we  suggest  that  the  observed  structure  In 
this  range  (1  -  10  km)  may  be  due  to  the  _E  x  B_  interchange  instability  due 
to  the  ambient  electric  fields,  or  neutral  winds,  or  a  combination  o' 
both.  We  have  not  considered  the  possible  origins  of  the  large  scale 
length  structures  (blobs,  patches,  etc.)  observed  at  high  latitudes  in  this 
work  (scales  sizes  >  100  kms ) ;  on  the  basis  af  this  work,  these  are 
produced  by  some  other  mechanism  than  the  interchange  process .  We  nave 
also  r.ot  attempted  to  in  elude  the  dynamics  of  these  large  structures  :r 


their  realistic  profiles.  The  interchange  process  discussed  here  level  ops 
over  a  time  scale  of  minutes  (-  102  sec)  whereas  the  long 
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variations  occur  over  a  period  of  hours,  thus  making  the  assumption  )'  i 
’static'  blob  valid  to  lowest  order. 
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In  conclusion,  we  have  argued  that  for  the  f inite-length  blobs  In  the 
high  latitude  F-region,  the  intermediate  scale  size  structure  responsible 
“or  satellite  scintillations  may  be  caused  by  the  _E  x  3.  interchange 
instability,  even  in  the  presence  of  a  conducting  E-region  (for  the  blob 
s^zes  of  -  300  km).  We  have  also  shown  that  “or  these  scale  sizes  (1  -  10 
km)  the  effect  of  ion- inertia  introduces  modest  reductions  in  the 
interchange  instability  growth  rates.  All  these  effects  (“inite  blob 
length,  ef'ect  of  inertia,  etc.)  however  act  to  reduce  the  growth  rate  of 
the  interchange  instability  at  long  scales  (>>  10  km)  severely,  thus  making 
it  unlikely  that  the  large  scale  sizes  would  be  locally  produced  by  the 
interchange  process  in  high-latitude  F-region. 
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APPENDIX  A:  LOCAL  THEORY 


We  Courier  analyze  the  perturbed  quantities  in  the  z-direction 
f(z)  -  f  exp  (ikzz)  for  the  local  approximation.  Then  (23)~(2^)  result  in 
the  dispersion  relation 


[y  +  ][y(1  +  -1-)  -  yJ  -  -  [YaDjkJ 


The  effects  of  finite  parallel  wavelength  are  included  on  the  right-hand 


side  in  the  above.  We  first  discuss  the  case  k? 
consider  the  case  when  k2  »  0. 


0,  and  later  would 


(a)  k. 


Setting  k  =»  0  in  (A1)  leads  to  the  dispersion  relation 


[y  +  k2d„][y(i  ♦  -2-)  -  yJ  -  o 


This  equation  has  solution  for  a  damped  hydromagnetic  mode,  Y  *  -  k^D^, 
[Chu  et  al.,  1978]  and  the  other  solution  refers  to  the  E  x  b  instability 
including  the  ion-inertial  effects  [Ossakow  et  al.,  1978] 

Y(l  ♦  -7-}  -  Y  -  0  (A3) 

vin  J 

The  two  limits  discussed  in  Ossakow  et  al.  [19"-]  are  the  collisional  and 
inertial  limits  for  the  solution  of  the  quadratic  in  (A3).  In  the 
collisional  limit,  v^n  >>  4yq,  one  obtains  the  well-known  result  of  local 
theory  growth  rate  of  the  E  *  3  instability  for  barium  clouds, 


In  the  inertial  limit,  the  growth  rate  of  the  E  *  3  instability  is  reduced 

from  the  collisional  case,  though  not  directly  proportional  to  a  decrease 

in  (as  might  be  expected  at  first),  but  proportional  to  /vin  instead 

[Ossakow  et  al . ,  1978].  Thus,  for  v.  <  4y„,  the  growth  rate  Is 

in  0 


Y  *  /Y0uln 


We  discuss  the  eollisional  and  inertial  limits  in  finite  k7 
approximation  for  the  electrostatic  and  electromagnetic  regimes. 

(i)  Electrostatic  regime  (k^D^  >>  y) 

In  an  electrostatic  approximation  the  dispersion  relation  (A1 ) 
reduces  to 


The  eollisional 
growth  rate  is 


,2 

,2  *  ^  (1  *<.!§)-  Vi„ 

y 


(A6) 


2  2 

limit  corresponds  to  4y.  <<  v.  (l  *  a  k  /k  1,  and  the 

0  in.  .  z  y‘ 


Y  -  Y  /(l  +  a  — ) 
k 

y 


(A7) 


Clearly,  the  parallel  dissipative  electron  motion  has  an  overall  effect  of 
reducing  the  growth  rate  in  this  limit. 

p  p 

The  inertial  limit  is  4Yn  >>  v.  fl  +  a  k  /k  )  and  the  growth  rate  in 

0  in  .  z  y 

this  case  is 


Y 


/Y  v, 
o  in 


(A8) 
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(ii)  Electromagnetic  regime  (y  >>  k^D^) 

In  this  case,  due  to  inductive  effects,  the  mode  couples  to  a 

* 

parallel  propagating  Alfven  wave  which  has  a  stabilizing  influence  on  the 
instability.  The  dispersion  relation  is  now 


v.  Y  -  ( Ynv,  -  V2k2) 

m  '  0  In  A  z‘ 


A9) 


The  growth  rate  in  the  eollisional  limit  is 


APPENDIX  3:  NONLOCAL  THEORY 


We  oar.  rewrite  the  nonlocal  dispersion  relation  (31)  Lr.  terns  of  tne 


normalized  parameters  as 


-  v . 

Y  -HI 


V  .  ,  V .  * 


tan.2k,  z. 


ir.  v/2 


-  Y  v .  /  v , 
m  i; 


ir.  ,1/2-  ir  ,*/ 


For  firite  size  alob?  and  at  higher  altitudes,  the  '’educed  growth 
rate  implies  Y  ■;  *  ar.d  ir.  this  limit,  ore  may  take  tan2k.  z-  *  2k,  z«.  Then 
the  dispersion  relatitr.  becomes 


Y[/a,  ♦  (l  +  — — )'/2  /a  ] 


*  ? ,  •  ,•-> 
k 

y 


[i  -  *  id  ♦  -!-)  *  /ri:  (i  *  H— ) 1 72 i ] 
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where  ap  -  ap  =■  ^./ve‘?Vp\  vin  *  Vtn/Y0’  In  the  ab0ve  we  have 

ignored  the  ion-inertial  effects  “or  the  bottomside  region  (1)  but 
retained  them  for  the  blob  and  the  topside  region  (2).  Now  we  write  down 
analytical  expressions  for  the  growth  rate  when  zQ  is  small. 


When  both  the  topside  and  bottomside  regions  and  the  blob  are 


collisicral ,  the  growth  rate  is 

^  A 

-  22  k 

Y  *  - Hi -  (33) 

Va,  +  /a,] 

where  a,  ar.d  a?  represent  the  “act  that  we  have  considered  different 
coil  is ichal ities  “or  all  the  three  regions.  If  ore  considers  the 
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situation  in  whioh  both  the  topside  and  bottomside  region?  may  be 
approximated  to  be  the  same,  then  one  obtains  the  result  presented  by 
Sperling  et  al.  [1981*] 

i 

i 

y  -  (z^kyj/Za  ' B4 ) 

l 

[  Case  II: 

» 

In  this  case,  the  bottomside  and  the  blob  are  assumed  coll is  tonal  but 
the  topside  is  considered  inertial  (v^n2  <  Yq],  then  the  growth  rate  is 
given  as 


2znk 

— 0  y 


(/al  +  /a2/Zin2] 


where  vln2  -  <  1. 


(35) 


Case  III: 

In  this  case,  the  bottomside  is  considered  colllslonal  but  the  ion- 
inertial  effects  are  included  for  both  the  topside  region  and  the  blob. 
Then  the  growth  rate  is 

y  -  (z0ky]/l(/ai  *  /  !2\  (36) 

vin2 


where 


